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Coupled  s-Wave  and  d-Wave  States  in  the  Heavy  Fermion 
Superconductor  ^Th^Be^ 
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239  Fronczak  Hall 

State  University  of  New  York  at  Buffalo 
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Abstract 

In  the  heavy- fermion  superconductor  ^h^Be^^ .  superconducting  states 
coexist  for  thorium  concentrations  0  <  x  <  0.06.  Assuming  s-wave  and  d-wave 
symmetries  for  these  states,  we  derive  a  Ginzburg -Landau  free  energy 
expression  which  couples  s-  and  d-wave  states  and  is  rotationally  invariant, 
in  contrast  to  the  free  energy  expression  proposed  by  Kumar  and  Wolfle  (Phys. 
Rev.  Lett.  59,  1954  (1987)).  We  discuss  in  detail  the  consequences  that 
follow  from  our  free  energy  relation.  In  particular,  we  predict  that  in  the 
above  system  there  are  two  eigenfrequencies  associated  with  the  dynamics  of 
phase  oscillations  (internal  Josephson  effect)  which  are  characteristic  of  the 
s-wave  and  d-wave  states. 


PACS  numbers:  74.70.Tx,  64.60.Cn,  74.20.De 


In  Che  heavy- fermion  system  ^Th^Be^.  different  superconducting  states 

have  been  found^  ^  to  coexist  for  impurity  concentration  x  in  the  range  0  <  x 

<  x^  (x^  =  0.06).  Two  superconducting  transitions  are  observed  in  specific 
12  3  4 

heat  ’  and  critical  field  ’  measurements.  For  undoped  UBe^.  the  higher 

temperature  transition  occurs  at  -  0.87  K,  whereas  the  lower  transition 

occurs  at  -  0.55  K.  Upon  doping  with  thorium  beyond  a  critical 

concentration  x^  =  0.017,  i.e.,  X2  >  x  >  x^,  the  high- temperature  and  low- 

temperature  superconducting  states  are  reversed  from  that  in  the  undoped 
3  4 

material.  ’  We  assume  the  two  superconducting  states  in  question  to  possess 
even-parity  spin-singlet  s-  and  d-wave  components.  We  are  aware  that  an 
interpretation  of  some  of  the  properties  of  the  superconducting  state  of  the 

Q 

UBe^j  system  in  terms  of  an  odd-parity  state  can  also  be  forwarded;  however, 
it  is  our  interest  to  study  the  relevance  of  s-  and  d-wave  coupling  in  this 
material . 

The  study  of  coupling  between  s-wave  and  d-wave  states  is  also  relevant 

7  9 

to  other  heavy-fermion  systems,  e.g.,  La-doped  UBe^  and  UPt^.  For  the 

high-Tc  superconductors,  the  relevance  of  s-wave  and  d-wave  coupling  is 

suggestive  from  the  microscopic  theories^ ’ ^  based  on  the  Hubbard  model. 

12 

High  resolution  X-ray  scattering  experiments  on  the  123  -  superconductor , 

¥8320^0^,  indicate  that  there  might  exist  a  special  interference  betwen  d- 

and  s-wave  states.  The  coupling  of  superconducting  states  is  reminiscent  of 

3 

the  classic  A-  to-B  transition  in  superfluid  He  which  involves  coupling 

between  different  spin-polarization  states  of  the  odd-parity  triplet  pairing 
r  •  1  13,14  rt  1  ,  »  ■  ,  ,  j  _  1 


of  the  quasi-particles. 


Superconductivity  due  to  d-wave  and  s-wave 


states  has  also  been  considered  by  Anderson  and  Morel  and  by  Mcrroin  and 


Stare . 


In  this  paper  we  develop  a  model  to  explain  the  important  electronic  and 


f?j(k)  -  Aj  (k)  exp(itfj)  ,  j  -  0,1,2 


(lb) 


magnetic  transitions  observed  in  the  xT^x^e^2  system.  The  higher 
temperature  transition  in  pure  UBe^  is  due  to  a  d-wave  state,  while  the  lower 
transition  is  due  to  an  s-wave  state  coexisting  with  the  d-wave  state.  We 
assume,  following  Kumar  and  Wolfle  (KW),^  that  thorium  doping  reduces  the 
transition  temperature  for  the  anisotropic  d-waves  due  to  impurity  scattering 
whereas  the  transition  temperature  for  the  isotropic  s -waves  is  hardly 
affected.  Thus,  for  x  >  x^,  the  d-wave  transition  temperature  T^)  is 

reduced  to  a  value  lower  than  that  of  the  s-wave  transition  temperature  Tq(- 


V- 


In  even-parity  spin-singlet  superconducting  states,  the  energy  gap 
function  A(k)  is  related  to  the  anomalous  thermal  average  <c£tc  of  the 
microscopic  theory,  c£t  being  the  electron  annihilation  operator  with  wave 
vector  1c  and  spin  t .  We  expand  A(lc)  as  a  linear  combination  of  the  basis 
function  set  {^.(k)},  where  <£0(k)  -  (Cx  +  cy  +  Cz)//3,  ^(k)  -  (Cy  -  Cx)  and 
“  (Cx  +  C  -  2Cz)//3  with  -  cos(k^a)  (i  -  x,y,z  and  a  the  lattice 
constant) : 


A(k) 


-  ^  00  (*> 


(la) 


j-o 


We  restrict  our  analysis  to  angular  momentum  states  £<  2  upon  which  the  basis  □ 
functions  (k)  can  be  replaced  by  the  spherical  harmonics  Yj (k) .  Then  j  -  0 
labels  the  s-states  and  j  -  1,2  label  the  two  d-state  components  of  the  E 
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irreducible  representation  of  the  octahedral  point  group,  0,  .  The  order 
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c«gA  q  +  AqAq  +  «2^2  +  ^2^2  +  "|,2^0^2  +  "^2^2^0^2COS  ^2  +  ”|r1^0^2cos^2 ’  (2a) 


parameter  amplitudes  A^ (k)  and  the  phases  6^  are  all  real.  Yg  is  a 

constant,  whereas  and  Y2  are  analogous  to  the  d  ^  2  anci  ^  2  or^^tals> 

x  -y  zL 

respectively,  of  atomic  physics.  The  symmetry  group  of  the  gap  function 
18  -  20 

is  G  x  R  x  U(l) ,  where  G  is  the  crystal  point  group,  R  is  the  time 

reversal  operator,  and  U(l)  is  a  gauge  transformation  group.  In  coordinate 

space,  the  Ginzburg- Landau  free  energy  of  the  superconducting  states  (relative 

to  the  normal  state)  can  be  written  as  F  -  J  d^r  (Fg  +  Fg) ,  where  the  free 

energy  density  F.  is  restricted  to  fourth  order  in  the  order  parameters  and  F„ 
L  G 

21 

contains  the  gradient  terms.  We  write,  Fg  -  F^  +  Fg,  where, 


and 


2  4  22  222  2  2 

Fg  -  a2^i  +  @2^1  +  "|,2^0^1  +  72^2^0^1COS  ^1  + 


+  7<SA^A2Cos2(0^-02)  '  7^4gA2A2  [2cos^2  +  cos(20^-0 ,) 


(2b) 


In  deriving  Eq.  (2)  we  have  adopted  the  notation  of  Kumar  and  Wolfle 

(KW)^  to  facilitate  comparison  with  their  work;  our  Eq.  (2a)  is  identical  to 

their  Eq.  (1).  Note  that  Eq.  (2a)  alone,  without  Eq.  (2b),  is  not 

rotationally  invariant,  since  the  d-wave  order  parameter  amplitudes  A^  and  A2 

transform  as  partner  components  of  the  irreducible  two-dimensional 

representation  E^  of  the  cubic  group  0^.  Equation  (2)  belongs  to  the  product 

3  2 

representation  Ag  x  E^.  Note  also  that  the  coupling  terms  AqA2  and  AqA^^ 
arise  from  the  fact  that  J  dll  Y^fc)  *  0  and  J  dfl  Yj  (k)Yg(fc)2  *  0, 
respectively;  these  two  terms  transform  into  each  other  under  symmetry 


.  ft*  A  V.WHaV  it 


operations  of  the  0^  group.  The  gauge  invariance  of  Eq .  (2)  follows  from  the 
fact  that  9^  and  9~  are  measured  relative  to  0q,  the  phase  angle  of  the  s-wave 
state . 

As  usual,  we  assume  that  a q  -  Aq(T-Tq)  and  a ^  "  Aj^-Tj),  where  Tq  and 
are  the  uncoupled  s-wave  and  d-wave  critical  temperatures,  respectively,  and 
that  the  other  coefficients  are  independent  of  temperature.  The  coefficients 
£q  and  ^2  are  assumed  positive  and  large,  relative  to  the  other  fourth-order 
coefficients,  to  insure  stability.  We  assume  the  condition  T2  >  Tq  for  the 
undoped  material,  whereas  the  reverse  is  true  for  the  doped  material  with  > 
x  >  x^  (X2  =*  0.06  and  x^  =<  0.017).  We  consider  the  following  seven  possible 
superconducting  states:  (I)  A^  0 ,  Aj  -  Aq  -  0,  (II)  *  0,  A^  -Aq  -  0, 
(III)  A1,  A2  *  0,  Aq  -0,  (IV)  Aq  **  0,  Ax  -  A2  -  0,  (V)  AQ ,  A±  *  0 ,  A2  -  0, 

(VI)  Aq  ,  A2  *  0,  A^  *■  0  and  finally  (VII)  AQ,  A^ ,  A2  *  0. 


For  a  transition  to  a  superconducting  d-state,  the  case  T2  >  TQ  and  T  « 
Tj  is  appropriate.  For  (a)  76  >  0  and  ^  2  “  ^  1  ~  state  I  or  state  III  is 

favored  depending  on  whether  7  >  2^  or  2^  -  |y|  >0.  respectively:  for  (b) 

76  <  0  and  9 2  “  8 y  ~  SCate  I  or  state  III  is  favored  depending  on  whether 
7(1+6)  >  2^2  or  2^2  '  | 7(1+6 ) |  >  0,  respectively.  Near  T2(>  Tq)  state  I  is 
stable  with  respect  to  non- vanishing  Aq  and  A^  since  these  two  order 
parameters  produce  a  second-order  effect  through  AqA2-  State  II  is  normally 
degenerate  with  state  I  in  the  absence  of  coupling;  however,  the  (s+d)- 
coupling  lifts  this  degeneracy  and  makes  state  II  unstable  with  respect  to 
state  VI.  With  a  decrease  of  temperature,  state  III  goes  over  to  state  VII 
since  coupling  to  Aq  becomes  important.  In  view  of  the  restrictions  on  the 
fourth-order  coefficients  imposed  in  the  previous  paragraph,  we  favor  state 
III  as  the  higher- temperature  d-state  for  the  undoped  material.  Taking  7  <  0, 


to 


r.!‘i 


1:1. 


7^>0,  72  <0,  6  <  0  and  <5 ^  <  0.  we  have  minimized  Eq.  (2)  numerically  in  the 
five -dimensional  space  of  the  order  parameters  and  the  relative  phase  angles. 
Our  plots  for  the  normalized  order  parameters  and  the  phase  angle  0^,  as  a 
function  of  temperature  are  given  in  Fig.  1;  9 ^  exhibits  a  constant  value  of 
rr/2. 

For  a  superconducting  transition  from  the  normal  to  the  s-wave  state,  as 

occurs  in  the  doped  material,  the  choice  Tq  >  with  T  =*  Tq  is  appropriate. 

Taking  all  other  parameters  exactly  the  same  as  in  Fig.  1,  state  IV  is 

favored.  With  a  decrease  of  temperature,  the  d-wave  states  get  populated  at  a 

temperature  T2  which  is  slightly  up-shifted  from  T2.  The  s-wave  and  d-wave 

order  parameters  coexist  giving  rise  to  state  VII.  Note  that  state  V  is 

unstable  with  respect  to  state  VII,  since  the  linear  term  involving  acts 

as  an  "applied  field"  and  lowers  the  free  energy.  Our  numerical  results 

appropriate  for  the  doped  materials  are  plotted  in  Fig.  2.  For  this  figure  we 

22 

find  that  state  VII  is  stable  with  respect  to  state  VI. 

The  d-wave  superconducting  states  have  some  characteristic  signatures 
such  as  power- law  temperature  dependencies  in  specific  heat,  muon  spin 
rotation,  nuclear  magnetic  resonance  and  ultrasonic  sound  attenuation 
experiments.  These  are  associated  with  vanishing  of  the  gap  function  along 
certain  symmetry  lines  or  planes  near  the  Fermi  surface.  Another 
characteristic  signature  is  the  dynamics  of  the  relative  phase  angles  0-^and  9^ 
due  to  (s+d) -coupling,  similar  to  the  internal  Josephson  effect  in  superfluid 
^He.^’^  The  particle  numbers  and  Nj  (relative  to  Nq)  are  conjugate  to 
the  phase  angles  0^and  9 ^  and  obey  the  commutation  relations  [  0  ^ ,  N  ^  ]  -  iS^ . 

We  generalize  Eqs.  (2)  to  include  the  chemical  potentials  and  linearize  the 
associated  Hamilton's  equations  about  0^  and  0^,  the  minimum  solutions.  The 
characteristic  phase  oscillation  frequencies  «  are  then  obtained  by  solving 


2  2 

the  secular  equation,  (a^-u>  Ka22'tt’  )  -  ai2a21  ”  w^ere  a- •  - 
2 

(3  i,j  -  1,2.  The  solution  of  the  characteristic  equation  is 

23 

analogous  to  that  of  a  double  pendulum  in  which  the  higher  frequency  normal 
mode  is  connected  with  a  symmetric  vibration  and  the  lower  frequency  one  with 
an  asymmetric  vibration.  For  solutions  corresponding  to  Fig.  1,  we  find  that 
the  high-frequency  mode  appears  at  the  first  transition,  whereas  the  low 
frequency  mode  appears  at  the  second  transition.  In  contrast,  both  the  normal 
mode  solutions  for  Fig.  2  appear  at  the  second  transition  (since  the  first 
transition  is  to  a  pure  s-wave  state,  there  is  no  coupling).  These 
charactertistic  phase  oscillations  thus  provide  a  method  to  identify  the 
superconducting  states.  At  low  temperatures,  the  high-frequency  fundamental 
vibration  can  couple  to  even- harmonics  (the  odd-harmonics  are  absent)  of  the 
low-frequency  fundamental  mode,  thus  producing  a  resonance.  Our  predictions 
in  this  regard  are  in  sharp  contrast  to  that  of  KW,  who  argued  for  a  single 
dynamic  frequency.  We  find  that  the  fundamental  frequencies  increase 
monotonically  as  T  -»  0,  and  we  do  not  find  any  evidence  for  non-monotonic 
behavior  as  suggested  by  KW. 

Finally,  we  would  like  to  comment  on  the  gradient  terms  F-,  of  the  free 

energy  density.  For  an  applied  magnetic  field  H  along  the  z-direction,  H  - 

(0,0,H),  we  choose  the  vector  potential  in  the  Landau  gauge,  A  -  (O.Ay,0),  Ay 

-  xH.  With  <&q  -  hc/2e,  we  obtain  the  invariant  terms  by  setting  up  covariant 

differential  operators.  For  the  sake  of  simplicity,  we  assume  a  one 

dimensional  coordinate  dependence  of  the  order  parameter  amplitudes:  rjj(r)  =* 

»jj(x)  ~  Aj(x)  exp(iSj),  j  -  0,1,2,  denote  derivatives  with  respect  to  x  by 
24 

primes  and  obtain  for  F^  -  G^y  +  G^ 


v'  'A  AA  A  . 


a0^o(A02  +  <2*/<VXA0]  +  |«2U2[A'z  +  A'z  +  (2  */#0)V(A^  +  a£)1  ,  (3 


*2*2,2,  J.  I™  It  TA>2  ,  a<2 


2  2  2  2 

‘r/>‘  a4' 


GE  “  X 


l^d{A^2-A[2+73  A^Aj  cos ( 0 2 -  0 1 )  +  (2^/*0)2A2[A2-A2-73A1A 


2  cos(»2-«1) k  + 


la2^sd[A0(^Al  cos*l  '  A2  cos(?2)  ‘  (2jt/$0) zA^A0(73A1  cos^  +  A2  cos«2)]  .  (3b) 


Here  $q  is  che  flux  quantum;  £q,  C2>^d  and  are  coherence  lengths.  We 

have  separated  the  terms  in  Eq.  (3)  in  the  same  spirit  as  that  of  Eq.  (2). 

Note  that  the  coefficients  of  the  £  .-terms  mix  the  s-  and  d-waves;  note  also 

sd 

2 

that  there  is  an  extra  d-wave  invariant  whose  coefficient  is  £  which  was  not 

considered  by  KW.  Employing  Eqs.  (2)  and  (3),  we  obtain  the  differential 

Ginzburg -Landau  (GL)  equations  and  use  them  to  derive  expressions  for  the 

upper  and  the  lower  critical  fields  H  and  H  ,  respectively. 

C2  C1 

The  upper  critical  field  is  the  highest  field  at  which  the  normal  state 

2 

free  energy  equals  that  of  the  mixed  state.  For  simplicity,  we  set  f  -  0  and 
obtain: 


H  -  (*0/2*)[c<c£  -  £*)  +  2?2£2]/[(c+2)(CA-^K2]  ,  (for  Fig.  1),  (4a) 

Hc  -  (*0/2,r)(1+2c')/t<^0+2c^2)2  *  12c'4]4  ’  (for  Fi6‘  2)  •  (4b) 


The  coefficients  c  and  c'  are  defined  as  Aq/A  and  A/Aq  respectively  (A  -  A^  - 

A2)  and  are  related  in  a  complicated  way  to  the  nonlinear  terms  of  F^.  Since 

2 

c  and  c'  are  less  than  one,  the  effective  coherence  length  is  dominated  by 
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or  as  the  case  may  be.  Setting  them  to  zero  we  recover  the  expected 

limits.  In  deriving  Eq.  (4)  we  have  found  the  ground  state  simple  harmonic 

oscillator  solutions  to  the  linearized  GL  equations  to  be  adequate. 

The  lower  critical  field,  i.e.,  the  field  at  which  the  Meissner  state  is 

thermodynamically  stable,  is  directly  proportional  to  the  super-electron 

density.  Consequently,  the  transition  to  the  mixed  (s+d)-wave  state  should  be 

accompanied  by  an  increase  in  the  temperature  coefficient  of  H  as  observea 

C1 

25 

in  Ref.  3.  By  assuming  the  large  ■*.(-  A/O  solution  for  H  and  the  London 

C1 

2  -  2 

form  for  the  supercurrent,  we  can  write  Hc  -  4>Qin  Ke£f/(2n^e£f) >  with  A  ^  = 

( Ag)  +  (A')'2  +  (A^)  * 2 ,  where  (Aq)  2  -  16*2c*0<;  2A2/$>2  ,  (A’)'2  -  327r2c*2£  2A2/4>2 


2S2  '  0 


and  (A^) ' L  -  (J3/2)  (i d/i^cos(ff  8 ^)  .  The  large  effect  of  the  relative 

_  2 

phase  difference  (8^~8y)  on  (A^)  may  explain  the  anomaly  in  *e££  near  the 

26 

lower  transition  temperature  reported  in  Ref.  3. 

In  summary,  we  have  derived,  from  general  symmetry  principles,  the 
Ginzburg- Landau  free  energy  for  the  mixed  x  E^  representation  of  the 
octahedral  point  group  relevant  for  pure  and  thorium-doped  UBe^-  By  detailed 
numerical  work,  we  have  obtained  the  nature  of  the  superconducting  states  that 
might  exist  in  this  system.  We  have  predicted  the  existence  of  a  symmetric 
and  an  asymmetric  normal-mode  phase  oscillation  frequency  characteristic  of 
the  coupled  states.  Finally,  we  have  examined  the  behavior  of  the  critical 
fields  for  this  system. 
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Figure  Captions 

1.  Temperature  dependence  of  order  parameters  (OP)  Aq  and  A  (-A^-^)  and 

' 

phase  angle  0^  for  critical  temperatures  T2/Tc  -  1  and  Tq/T^  ”  0.75  and 
j  phase  angle  6^  -  ir/2.  This  figure  is  appropriate  for  undoped  UBe^-  The 

x-axis  and  the  left-hand  side  y-axis  variables  are  in  dimensionless  units. 


2.  Same  as  in  Figure  1,  but  with  Tq/Tc  -  1  and  T2/TQ  -  0.75.  This  figure  is 
appropriate  for  U^^Th^Be^^- 
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